What Is Claimed Is; 

1 1 . A method for using a computer system to solve an unconstrained 

2 interval global optimization problem specified by a function/ wherein/is a scalai 

3 function of a vector x = (x h x 2 , x 3 , ... x n ), the method comprising: 

4 receiving a representation of the function/at the computer system; 

5 storing the representation in a memory within the computer system; and 

6 performing an interval global optimization process to compute guaranteed 

7 bounds on a globally minimum value of the function fix); 

8 wherein performing the interval global optimization process involves, 

9 applying term consistency over a subbox X, and 

! o excluding any portion of the subbox X that violates term 

1 1 consistency. 

1 2. The method of claim 1 , wherein applying term consistency 

2 involves: 

3 symbolically manipulating an equation within the computer system to 

4 solve for a first term, g(x'), thereby producing a modified equation g(x') = h(x), 

5 wherein the first term g(x') can be analytically inverted to produce an inverse 

6 function g'(y); 

7 substituting the subbox X into the modified equation to produce the 

8 equation g(X') = h(X); 

9 solving for X ' = g '\h(X)); and 

I o intersecting X ' with the subbox X to produce a new subbox X + ; 

I I wherein the new subbox X + contains all solutions of the equation within 

12 the subbox X, and wherein the size of the new subbox X + is less than or equal tc 

1 3 the size of the subbox X: 
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3. The method of claim 1, wherein performing the interval global 
optimization process involves: 

keeping track of a least upper bound f_bar of the function f(\); 
removing from consideration any subbox for which f(x) >f_bar\ 
applying term consistency to the inequality f{\) <f_bar over the subbox X; 

and 

excluding any portion of the subbox X that violates the inequality. 

4. The method of claim 1, wherein performing the interval global 
optimization process involves: 

determining a gradient g(x) of the function /x), wherein g(x) includes 
components g,(x) (/=1, ...,«); 

removing from consideration any subbox for which g(x) is bounded away 
from zero, thereby indicating that the subbox does not include a global minimum 
of/[x); and 

applying term consistency to each component g/(x)=0 (z'=7, ...,n) of g(x)=0 
over the subbox X; and 

excluding any portion of the subbox X that violates a component. 

5. The method of claim 1, wherein performing the interval global 
optimization process involves: 

determining diagonal elements //,,(x) (/=1, ... t ri) of the Hessian of the 
function/(x); 

removing from consideration any subbox for which a diagonal element of 
the Hessian is always negative, which indicates that the /is not convex and 
consequently does not contain a global minimum within the subbox; 
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applying term consistency to each inequality #„(x) z 0 (*'=1, ...,«) over the 
subbox X; and 

excluding any portion of the subbox X that violates an inequality. 

6. The method of claim 1 , wherein performing the interval global 
optimization process involves performing the Newton method, wherein 
performing the Newton method involves: 

computing the Jacobian J(x,X) of the function / evaluated as a function of 
a point x over the subbox X; 

computing an approximate inverse B of the center of J(x,X); and 

using the approximate inverse B to analytically determine the system 
Bg(x) 5 wherein g(x) is the gradient of the function/x), and wherein g(x) includes 
components g,(x) (/=/, ...,«); and 

applying term consistency to each component (Bg(x)), = 0 (/=1, ...,ri) for 
each variable x, (i=l, ...,ri) over the subbox X; and 

excluding any portion of the subbox X that violates a component. 

7. The method of claim 1 , further comprising terminating attempts to 
further reduce the subbox X when: 

the width of X is less than a first threshold value; and 

the magnitude offlX) is less than a second threshold value. 

8. A computer-readable storage medium storing instructions that 
when executed by a computer cause the computer to perform a method for using a 
computer system to solve an unconstrained interval global optimization problem 
specified by a function/ wherein /is a scalar function of a vector x = (xj } %2, x$, ... 
x n ), the method comprising: 
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receiving a representation of the function / at the computer system; 
storing the representation in a memory within the computer system; and 
performing an interval global optimization process to compute guaranteed 
bounds on a globally minimum value of the function/(x); 

wherein performing the interval global optimization process involves, 
applying term consistency over a subbox X, and 
excluding any portion of the subbox X that violates term 
consistency. 

9. The computer-readable storage medium of claim 8, wherein 
applying term consistency involves: 

symbolically manipulating an equation within the computer system to 
solve for a first term, g(x'), thereby producing a modified equation g(x') = h(x), 
wherein the first term g(x') can be analytically inverted to produce an inverse 
function g\y); 

substituting the subbox X into the modified equation to produce the 
equation g(X') = /z(X); 

solving for X J = g~\h(X)); and 

intersecting X ' with the subbox X to produce a new subbox X + ; 

wherein the new subbox X + contains all solutions of the equation within 
the subbox X, and wherein the size of the new subbox X + is less than or equal to 
the size of the subbox X. 

10. The computer-readable storage medium of claim 8, wherein 
performing the interval global optimization process involves: 

keeping track of a least upper bound fjbar of the function/(x); 
removing from consideration any subbox for which/(x) >f_bar\ 

27 

Attorney Docket No. SUN-P7076-SPL Inventors: Walster et al. 

ARPCAMY DOCUMENTS\SUN MICROSYSTEMS\SUN-P7076-SPL\SUN-P7076 APPLICATI0N.DOC 




applying term consistency to the inequality j(x) <f_bar over the subbox X; 

and 

excluding any portion of the subbox X that violates the inequality. 

1 1 . The computer-readable storage medium of claim 8, wherein 
performing the interval global optimization process involves: 

determining a gradient g(x) of the function^*), wherein g(x) includes 
components g,(x) (7=1, .. .,/?); 

removing from consideration any subbox for which g(x) is bounded away 
from zero, thereby indicating that the subbox does not include a global minimum 
of/00; an <3 

applying term consistency to each component g{x)=0 (7=7, ...,«) of g(x)=0 
over the subbox X; and 

excluding any portion of the subbox X that violates a component. 

12. The computer-readable storage medium of claim 8, wherein 
performing the interval global optimization process involves: 

determining diagonal elements //,,(x) (7=1, ...,«) of the Hessian of the 
function/(x); 

removing from consideration any subbox for which a diagonal element of 
the Hessian is always negative, which indicates that the /is not convex and 
consequently does not contain a global minimum within the subbox; 

applying term consistency to each inequality //,,(x) ^ 0 (7=1, ...,/?) over the 
subbox X; and 

excluding any portion of the subbox X that violates an inequality. 
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13. The computer-readable storage medium of claim 8, wherein 
performing the interval global optimization process involves performing the 
Newton method, wherein performing the Newton method involves: 

computing the Jacobian J(x,X) of the function / evaluated as a function of 
a point x over the subbox X; 

computing an approximate inverse B of the center of J(x,X); and 

using the approximate inverse B to analytically determine the system 
Bg(x), wherein g(x) is the gradient of the function /(x), and wherein g(x) includes 
components g,(x) (i=l, ...,/?); and 

applying term consistency to each component (Bg(x))/ = 0 (/=1, ...,») for 
each variable x, (/=1, ...,«) over the subbox X; and 

excluding any portion of the subbox X that violates a component. 

14. The computer-readable storage medium of claim 8, wherein the 
method further comprises terminating attempts to further reduce the subbox X 
when: 

the width of X is less than a first threshold value; and 

the magnitude of XX) is less than a second threshold value. 

15. An apparatus that solves an unconstrained interval global 
optimization problem specified by a function / wherein /is a scalar function of a 
vector x = (xj, X2, X3, ... x„), the apparatus comprising: 

a receiving mechanism that is configured to receive a representation of the 
function /; 

a memory for storing the representation; and 
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an interval global optimization mechanism that is configured to perform 
an interval global optimization process to compute guaranteed bounds on a 
globally minimum value of the function /(x); 

a term consistency mechanism within the interval global optimization 
mechanism that is configured to, 



16. The apparatus of claim 15, wherein the term consistency 
mechanism includes: 

a symbolic manipulation mechanism that is configured to symbolically 
manipulate an equation within the computer system to solve for a first term, g(x'), 
thereby producing a modified equation g(x') = h(x), wherein the first term g(x') 
can be analytically inverted to produce an inverse function g'\y); 

a solving mechanism that is configured to, 



an intersecting mechanism that is configured to intersect X' with the 
subbox X to produce a new subbox X + ; 

wherein the new subbox X + contains all solutions of the equation within 
the subbox X, and wherein the size of the new subbox X + is less than or equal to 
the size of the subbox X. 

1 7. The apparatus of claim 1 5, 

wherein the interval global optimization mechanism is configured to, 



apply term consistency over a subbox X, and to 
exclude any portion of the subbox X that violates term 



consistency. 



substitute the subbox X into the modified equation to 
produce the equation g(X') = /z(X), and to 
solve for X' = gl '\h{X))\ and 
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keep track of a least upper bound f_bar of the function/x), 

and to 

remove from consideration any subbox for which 



wherein the term consistency mechanism is configured to, 

apply term consistency to the inequality fix) <,f_bar over 

the subbox X, and to 

exclude any portion of the subbox X that violates the 

inequality. 

1 8 . The apparatus of claim 1 5 , 

wherein the interval global optimization mechanism is configured to, 

determine a gradient g(x) of the fimctiony(x), wherein g(x) 
includes components g,(x) (/=1, ...,/?), and to 

remove from consideration any subbox for which g(x) is 
bounded away from zero, thereby indicating that the subbox does 
not include a global minimum of/(x); and 
wherein the term consistency mechanism is configured to, 

apply term consistency to each component gj(x)=0 
(/=7, ...,ri) of g(x)=0 over the subbox X, and to 

exclude any portion of the subbox X that violates a 
component. 

1 9. The apparatus of claim 1 5, 

wherein the interval global optimization mechanism is configured to, 
determine diagonal elements H u (x) (z=l, ...,n) of the 
Hessian of the function fix), and to 
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Ax) >f_bar- 
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remove from consideration any subbox for which a 
diagonal element of the Hessian is always negative, which 
indicates that the /is not convex and consequently does not contain 
a global minimum within the subbox; and 
wherein the term consistency mechanism is configured to, 

apply term consistency to each inequality ////(x) ^ 0 
(/=1, ...,«) over the subbox X 5 and to 

exclude any portion of the subbox X that violates an 
inequality. 

20. The apparatus of claim 15, further comprising a Newton 
mechanism within the interval global optimization mechanism that is configured 
to: 

compute the Jacobian J(x,X) of the function / evaluated as a function of a 
point x over the subbox X; 

compute an approximate inverse B of the center of J(x,X); and to 

using the approximate inverse B to analytically determine the system 
Bg(x), wherein g(x) is the gradient of the function/x), and wherein g(x) includes 
components g,(x) (z=l, ...,n); 

wherein the term consistency mechanism is configured to, 



apply term consistency to each component (Bg (x)), = 0 
(/=/, ...,«) for each variable x, (z'=7, ...,«) over the subbox X, and to 

exclude any portion of the subbox X that violates term 
consistency. 
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21 . The apparatus of claim 15, further comprising a termination 
mechanism that is configured to terminate attempts to further reduce the subbox X 
when: 

the width of X is less than a first threshold value; and 

the magnitude of/(X) is less than a second threshold value. 
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